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ABSTRACT 

Enigmatic transitions between spin-up and spin-down have been observed in several 
X-ray pulsars accreting matter via an accretion disk. In these transitions, the torque 
changes sign but remains at nearly the same magnitude. It has been noted previously 
that alternating prograde and retrograde disk flows would explain many features of 
the torque reversals, although it has been unclear how a stable retrograde disk could 
be formed. We suggest that the reversals may be related to the disk at times being 
warped to such an extent that the inner region becomes tilted by more than 90 degrees. 
This region would thus become retrograde, leading to a negative torque. Accretion 
disk models can show such behavior, if account is taken of a warping instability due to 
irradiation. The resulting 'flip-overs' of the inner parts of the disk can reproduce most 
characteristics of the observations, although it remains unclear what sets the timescale 
on which the phenomenon occurs. If this model were correct, it would have a number 
of ramifications, for instance that in the spin-down state the X-ray source would mostly 
be observed through the accretion disk. 
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1. Introduction 

Long-term, continuous monitoring by the BATSE 
all-sky monitor on the Compton Gamma Ray Obser- 
vatory has led to a qualitative change in our picture 
of the spin-frequency behavior of accre ting X-ray pul- 
sars (Bildsten et al. 1997, hereafter |b"97| ). Of the 
four well-measured persistent sources thought to ac- 
crete by way of an accretion disk, all display sudden 
transitions between episodes of steady spin-up and 
spin-down. The spin change rate (i.e., the absolute 
value of the pulse frequency derivative) is nearly equal 
for spin-up and spin-down, and it is comparable to 
the torque expected if all the angular momentum of 
the accreting gas is deposited at the magnetospheric 
boundary and transferred to the neutron star. In 
at least some of these systems, however, Roche-lobe 
overflow is thought to occur, and hence the matter in- 
serted into the accretion disk has only one sense of an- 
gular momentum; symmetric torque reversals would 
thus not be expected. 

Previous accretion torque models, in which the net 
torque could be positive or negative depending on the 
mass accretion rate M (e.g., phosh fc Lamb 1979 ) 



cannot easily reproduce the observations (see Nelson 



t al. 1997). We only list the main arguments here: (i) 



one needs step- wise changes in M to produce distinct 
spin-up and spin-down states - this seems unlikely; 
(ii) one would expect changes in M to be reflected 
in the X-ray luminosity Lx - there are variations in 
ix, but these appear uncorrelated with spin- up/down 
state; and (iii) one expects at all times a positive cor- 
relation between torque and luminosity - for GX 1+4 
an anti-correlation is observed during its spin-down 
state ( |Chakrabarty et al. 19971] ). 

Suggestions about a possible cause for the torque 
revers als have been made by Yi , Wheeler & Vishniac 
( [1997D and Nelson et al. ( [19971) . Yi et al. suggested 
that the reversals were due to small changes in M 
around a critical value at which the system changes 
from a primarily Keplerian flow to a substantially 
sub-Keplerian, radially advective flow. This addresses 
points (i) and (ii), but not point (iii). 

Nelson et al. ( |1997 ) explored the possibility of hav- 
ing systems in which nothing changes except the sense 
of rotation of the disks. They found that this would 
explain the observations very well. They quoted 
a suggestion by Makishima et al. (1988), that in 
GX 1+4 one might have accretion from a wind and 
thus form a retrograde disk more easily. As Nelson et 



al. noted, however, it is very hard to imagine how a 
stable retrograde disk could form, especially in the ul- 
tracompact binary 4U 1626—67, for which all indica- 
tions are that mass transfer is by Roche lobe overflow 
from a very low-mass companion. 

Here, we suggest a modification of the Nelson et 
al. picture, viz., that it is only the inner part of the 
disk that is changing its sense of rotation, as a con- 
sequence of very strong warping of the disk. Accre- 
tion disks are unstable to warping if lit strongly by a 
central radiation source (Pringle 1996). In numerical 
simulations which include such irradiation, the disk 
can sometimes become more and more warped, un- 
til the inner part has become inclined by well over 
90 degrees (Pringle 1997; see Fig. 6 of that paper 



for illustrations of a disk with its inner parts flipped 
over). 

In Section ^, we give an updated summary of the 
observations. In Section ||, we briefly discuss the 
warping process and describe simulations done specif- 
ically for X-ray binaries. We proceed to make a qual- 
itative comparison with the observations, and to dis- 
cuss ramifications. We summarize our conclusions in 
Section [|. 

2. Summary of the Observations 

BATSE has monitored five persistent, disk-fed ac- 



creting X-ray pulsars for nearly 7 years (B97). Here, 
we summarize the observations for four of these. We 
exclude Her X-l, since it cannot be monitored con- 
tinuously, but only during the high states of its 35-d 
precession cycle. 

4U 1626—67 is a 7.66s pulsar in an ultracompact 
42min binary. It has Lx — 10 37 ergs _1 (Chakrabarty 



1998). X-ray timing measurements suggest we are 
viewing the binary nearly face-on. The < 0.1 Mq 
mass donor probably fills its Roche lobe. The torque 
history is reviewed in detail by Chakrabarty et al. 
(1997a). The pulsar underwent steady spin-up at 
a mean rate of 0.85 pHzs -1 during 1977-1990, and 
steady spin-down at a mean rate of — 0.72pHzs _1 
since that time. The 1990 reversal itself was not ob- 
served, but evidently occurred on a time scale <J 1 
month. Both the spin-up and the spin-down had a 
significant quadratic trend. 

The X-ray luminosity did not change abruptly at 
the torque reversal, although there is evidence for 
a gradual, factor of 2 decay over an interval of 15 
years. However, the long-term (years) 1-20 keV X- 
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ray spectral shape did change, from dN/dE oc E 15 
during spin-up to dN/dE oc E-° A during spin-down 



(Ve ughan & Kitamoto 1997), while the > 20keV spec 



trum remained unchanged. During spin-up, the pul- 
sar showed strong flares on ~ 1000 s time scales, but 
these flares are absent or much less frequent during 
spin-down (Chakrabarty et al. f998, in preparation). 

Quasi-periodic oscillations (QPO) were observed 
bo th during spin-up fO.041 Hz; Bhinoda et al. 199C \ 



an d spin-down (0.048 Hz; Angclini et al. 1995). The 



detailed sideband structure of the X-ray pulsations 
and the QPO during spin-down suggests that there 
is X-ray reprocessing or scattering in material on 
a prograde orbit, likely at small radius (Kommers, 
Chakrabarty, & Lewin 1998). 

GX 1+4 is a 2min pulsar in a wide (P rb <; 100 d) 
binary with an M giant, and its accretion disk is prob- 
ably fed by the M giant's wind (Chakrabarty fc Roche 



19S0). Most likely, L x ^ 10 37 ergs -1 , although a 
larger distance corresponding to Lx — 10 38 ergs _1 
cannot be not ruled out. The torque and flux historie s 
are reviewed in detail by Chakrabarty et al. ( 1997b ). 
The source was bright throughout the 1970s and spin- 
ning up rapidly at a mean rate of +6.0 pHz s _1 . EX- 
OS AT observations in 1983 and 1984 failed to detect 
it. When detected again at a low luminosity in 1987, 
it was spinning down at a mean rate of — 3.7pHz s _1 
(1989-1991). BATSE has observed continued spin- 
down since 1991, apart from a short spin-up episode 
and an episode in which the source was not detectable. 
The torque reversals observed by BATSE were re- 
solved in 5-day-averaged spin measurements. 

The long-term mean X-ray luminosity was brighter 
during the 1970s spin-up interval (and also during 
the BATSE spin-up event) than during the spin-down 
intervals in the 1980s and 1990s. However, during 
the extended spin-down interval observed by BATSE, 
an anti-correlation between torque and X-ray lumi- 
nosity was found: increased spin-down for increased 
X-ray flux. The luminosity during spin-down occa- 
sionally flared to values similar to those measured 
during spin-up in the 1970s. The X-ray spectrum 
did not change drastically between the 1970s and the 
late 1980s ( Sakao et al. 1990| ), although considerable 



variability in the hydrogen column density has been 
observed. 

Cen X-3 is a 4.8 s pulsar in a 2.1 d binary with an 
OB supergiant. It has Lx — 10 37 ergs _1 . Cen X-3 
has shown a long-term spin-up trend of 1.6pHzs _1 
since 1972, but BATSE observations resolve this trend 



into alternating episodes of spin-up and spin-down on 
10-100 day time scales ( |Finger, Wilson fc Fishmaii 
1994). The torque distribution is bimodal, with pre- 
ferred spin- up/spin-down rates of +7 and — 3pHz s _1 
( B97| ). The torque reversals are unresolved in the 
BATSE measurements and must occur in < 10 d. 
Over a factor of 6 range in luminosity, no correlation 
between torque and luminosity is found, and the mean 



flux is comparable in the two torque states (Vaughan 
et al. 1998). The spectral shape, pulse profile, and 



pulsed fraction arc correlated with luminosity. 

O40 1657-415 is a 38s pulsar in a 10.4 d bi- 
nary with an unidentified companion which is inferred 
to be an OB supergiant. It has Lx — 10 37 ergs _1 



(Chakrabarty et al. 1993). There is no direct evi- 
dence for an accretion disk, but the timing properties 
of the pulsar strongly suggest disk accretion. Like 
Cen X-3, the pulsar has shown a long-term spin-up 
trend since 1979 of ~ 1 pHz s _1 , but BATSE observa- 
tions show alternating episodes of spin-up and spin- 
down on 10-100 d tim e scales, with typical rates of +7 
and —2 pHz s _1 ( B97 ). The torque reversals are unre- 
solved in the BATSE measurements. No correlation 
between torque and luminosity is observed, although 
this has only been tested over a small range (factor of 
2) in luminosity flKoh 1998| ). 



3. Warped disks 

Wijers and Pringle (1997; hereafter WP97) have 
investigated in detail the behavior of irradiated accre- 
tion disks in binary stars subject to internal viscous 
forces, radiation reaction due to illumination from the 
central source, and the tide from the mass donor. The 
main aim was to see whether warped disks could form, 
and whether their precession could lead to the super- 
orbital periodicities observed in many X-ray binaries. 
Here, we summarize the results relevant to the torque 
reversals, and compare them to the observations. 

Whether the disk becomes unstable to warping de- 
pends on the parameter F*, given by (Eq. 18, WP97) 



F, = 0.127 t] ~~~ ^ 



.4 



1.4 M, 
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0.1 0.01 V 10 n cm Af NS 

where r\ is the ratio of the (R, z) and (R, 4>) vis- 
cosities (see Pringle 1996 ), e = Lx/Mc 2 the accre- 
tion efficiency, A the disk's X-ray albedo, Mns the 
neutron-star mass, and i? ou t the outer disk radius. 
WP97 find that for L* <; 0.1, the disks are stable; for 
0.1 <; F t <; 0.15, the outer parts become warped and 
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precess retrogradely, but the inner parts remain flat; 
for 0.f5 <; F * <; 0.2, the inner parts become warped 
as well, precessing progradely; and for F* > 0.2, the 
disks start to behave somewhat chaotically, precess- 
ing less steadily, and often having inner parts tilted 
by more than 90 degrees. 

The numerical value of F* in a real system is quite 
uncertain, mostly because the disk albedo is poorly 
known (see WP97), but for typical X-ray binary pa- 
rameters, one will be close to the unstable regime. We 
will assume that F» varies, and associate the switch 
between spin-up and spin-down state with a switch 
between a prograde and retrograde inner disk (for 
prograde pulsar rotation). What could drive the vari- 
ations is unclear, as F* does not depend directly on 
quantities such as M which are usually thought to 
vary. 

The one parameter which might lead to a switch 
between two states is the disk albedo. As discussed by 
WP97, the albedo likely depends on the angle of inci- 
dence, with (1 — A) < 0.1 for grazing and (1— A) ~ 0.4 
for normal incidence. One could imagine a situation 
in which a disk is stable, until a random excursion 
changes the tilt such that F* is raised above the criti- 
cal value, and the disk starts to become warped. This 
would increase (1 — A) further, leading to further 
growth, etc. It remains problematic, however, how 
one would reset the system. 

3.1. Inner-disk tilt 

An example of the time evolution of the disk tilt 
angle for F* = 0.2 is given in Fig. |l|. It shows that only 
the inner disk becomes tilted by more than 90 degrees 
(i.e., retrograde). Reversing the sense of rotation in 
the inner disk would lead to torque reversals similar 
to those observed, just as reversing the whole disk 
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would in the picture of Nelson et al. (1997). A retro- 
grade inner disk could also lead to the anti-correlation 
between torque and luminosity observed in the spin- 
down state of GX 1+4. 

In the model in Fig. [l], as in most of the simula- 
tions we tried, the average inner disk tilt does not 
reach 180°. In contrast, for most values of F* leading 
to prograde inner disks, the tilts are close to 0°. If 
generally true, then for a roughly spherical magne- 
tosphere the spin-change rate in the spin-down state 
should be systematically lower than that in the spin- 
up stateQ. For all four sources, the spin-down rate is 
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Fig. 1. — Disk tilt as a function of time at var- 
ious radii. The top panel shows the full simula- 
tion, and the bottom panel a part of it on an ex- 
panded scale. For the simulation, we used F* = 0.2, 
and matter enters at two thirds the outer radius 
(see Wijers & Pringlc 1997). In units of the outer 



radius, the radii shown are at 0.0083, 0.083, 0.47, 
and 0.87 (solid, long-dash, short-dash, dotted, respec- 
tively). For 4U 1626-67, Cen X-3, OAO 1657-415, 
and GX 1+4, the outer radii are approximately 0.18, 
2, 5, and 50 x 10 11 cm, respectively. For M = 
lO _9 M0yr _1 , one unit of time equals 0.3, 5, 17, and 
300 d, respectively. 



indeed smaller, by a factor 1.2 to 4. 

Variations around the mean of about 30° (rms) 
in the inner disk tilt angle occur with a character- 
istic period of ~ 4 time units. This corresponds to 
about 10 5 s in 4U 1626-67, 10 6 s in Cen X-3, 10 7 s in 
OAO 1657-415, and 10 8 s in GX 1+4. (Note, how- 
ever, that F* = 0.2 is not necessarily appropriate 
for all systems.) These variations result in torque 
fluctuations, which should contribute to the torque 



tation axis due to the non-parallel component of the inner-disk 
angular momentum are likely small. The effects are reduced by 
precession for GX 1+4 and 4U 1626—67, and by randomness 
in inner disk orientation for Cen X-3 and OAO 1657—415. 



*We note that changes in the orientation of the neutron-star ro- 
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noise observed in X-ray pulsars^. For Cen X-3 and 
OAO 1657—415, this may not be very relevant, as 
the characteristic period is comparable to the dura- 
tion of a spin-up or spin-down state. However, for 
4U 1626-67 and GX 1+4, it is much shorter than the 
state duration, and we can compare the noise induced 
by the fluctuations with the observations. 

For GX 1+4, the tilt variations correspond to a 
torque noise power of ~ 10~ 16 Hz 2 s~ 2 Hz -1 at 10~ 8 Hz. 
This is similar to what is observed (B97). For 4U 1626—67, 
we need to integrate the tilt variations over 10 7 s, the 
shortest timescale on which the torque noise power 
has been measured. We find 
somewhat larger than the observed < 
at 10 _7 Hz (Chakrabarty et al. 1997a), but perhaps 



not inconsistent given the uncertainties in what a real 
disk would do. 

3.2. Observing through the disk 

An interesting consequence of associating spin- 
down with a flipped-over state is that one should be 
viewing the neutron star through the disk at least 
part of the time. We estimate the surface density us- 
ing a steady, unwarped Shakura-Sunyaev disk. For 
all four systems, M ~ 1(T 9 M© yr" 1 . The typi- 
cal radius where the line of sight crosses the disk 
is where the disk tilt is 90°, about 0.25i? O ut in our 
simulations. In physical units, this ranges between 
- 10 9 - 5 and - 10 12 cm (in 4U 1626-67 and GX 1+4, 
respectively). With Mjjs = XAMq and viscosity 
parameter a, we find a surface density of 50a -4 / 5 
and 0.7a -4 / 5 g cm" 2 , respectively (Eq. 5.45 of Frank 



King, fc Rainc 1992| ). This corresponds to column 



densities 7V H of about 10 25 a~ 4/5 and 10 23 a- 4/5 cm"', 
respectively. The disk central temperature due to vis- 
cous heating is 10 4 5 a~ 1 / 5 and 10 3 a -1 / 5 K, respec- 
tively, and the ionization parameter £ = Lx/nR 2 a 
few a 7 / 10 (cgs) for both cases. Hence, little ioniza- 
tion is expected, and the disk would be opaque to 
low-energy X-rays for the dense case. 

In a twisted disk, the extra viscosity associated 
with the misalignment between neighboring annuli 
will lower the column density somewhat. Compared 
to a similar disk with no irradiation, the column den- 
sity in the warped disk shown in Figure ^ is lower by 
a factor 5 in the outer parts of the disk, and by a 
factor 2.5 in the inner parts. 



5 The two contributions usually considered are M vari ations and 
internal torques; see Lamb, Pines, & Shaham (|197q). 



Good data on possible spectral changes associated 
with the torque reversals exist only for 4U 1626—67 
(Sect. |2j). The observed decrease and hardening of the 
1-20 keV flux in that source may reflect a line of sight 
through the disk. The spectrum changed mostly be- 
low ~ lOkeV, suggesting that absorption might play 
a role. Vaughan & Kitamoto (1997), however, found 
that the spectrum could not be fit well by simply in- 
creasing the column density Ah; for an Ah sufficient 
to explain the reduction at close to 10 keV, one would 
expect no soft X-rays whatsoever. 

Increasing Ah may not be quite appropriate, as (i) 
matter in the disk is likely neither cold nor homoge- 
neous; (ii) scattering might be important, as the disk 
subtends a large solid angle; and (iii) the obscuration 
will vary will time, due to precession of the inner disk 
- for suitable inclinations, the source may even be 
unobscured part of the time. Of course, given free- 
dom to change all these, one likely can fit almost any- 
thing. Vaughan & Kitamoto could fit the 4U 1626—67 
spectrum assuming partial covering, where 86% of the 
source is obscured with Ah = 10 23 8 cm -2 and the re- 
mainder unobscured. If this were the correct model, it 
might reflect either spatial or temporal covering (the 
spectrum was a long-term average). 

3.3. Problems 

There are several difficulties in applying the warped 
disk model to the torque reversal observations in de- 
tail. The primary one concerns the time scales be- 
tween the reversals. In those simulations where a 
retrograde inner disk arose, its mean inclination re- 
mained relatively stable (and > 90°) after its initial 
growth, with only short excursions to prograde flow. 
This is at odds with the observations, which show 
roughly equal times in spin-up and spin-down. 

Furthermore, each of the four systems considered 
has an aspect that is not easily explained in the retro- 
grade disk picture. In 4U 1626—67, there is QPO evi- 
dence for prograde motion during spin-down, likely at 
small radius (Sect. ^J). In GX 1+4, the torque rever- 
sals observed by BATSE were clearly accompanied by 
luminosity changes. Finally, in both high-mass X-ray 
binaries (Cen X-3 and OAO 1657-415), the lack of a 
clear correlation between spin-change rate and lumi- 
nosity (in both spin-up and spin-down states) seems 
to confound all the existing models. 

A number of potentially important effects have 
been neglected in the WP97 simulations. These in- 
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elude: the dependence of the albedo on angle of in- 
cidence (and the resulting variations of F*), the in- 
creased disk scale height and thus reduced density 
due to irradiation (resulting in larger £ and possibly 
different disk opacity), and the varying "splash point" 
where the accretion stream - which is in the orbital 
plane - hits the precessing disk (resulting, e.g., in 
inward advection in addition to viscosity). Account- 
ing for these properly may considerably change the 
detailed predictions for warped disks. However, the 
existence of warps and transitions to retrograde flow 
under certain conditions seem a general feature. 

4. Conclusions 

It has been previously suggested that retrograde 
disk flows would be an elegant explanation for torque 
reversals in disk- fed X-ray pulsars. It has been an 
open question, however, how such a flow would be 
created when mass is transferred by Roche-lobe over- 
flow and thus enters the accretion disk in a prograde 
orbit. We have pointed out that strong disk warping 
may produce a retrograde flow close to an accreting 
neutron star from an initially prograde flow in the 
outer disk. 

While this simple picture allows us to understand 
some of the characteristics of the observations, it does 
not match all of them. This could be because the 
models currently neglects many potentially important 
effects. Our main point here, however, is that quasi- 
stable retrograde disk flows are physically plausible, 
and that they may thus be (part of) what causes the 
spin-down states in persistent disk-fed X-ray pulsars. 

Despite the uncertainties in the models, it remains 
true generally that if the spin-down state is due to a 
flipped-over inner disk, it should be associated with 
enhanced absorption. One would expect relatively 
strong absorption edges, especially of iron, which 
could be looked for with high-resolution spectroscopy. 
The absorption should vary with time as the inner 
disk precesses (although the precession is not neces- 
sarily very coherent). With more disk surface exposed 
to the central X-ray source, we also would expect in- 
creased scattering and stronger fluorescence lines. 
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